Influence of Strain on the Kinetics of Phase Transitions in Solids 
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We consider a sharp interface kinetic model of phase transitions accompanied by elastic strain, 
together with its phase-field realization. Quantitative results for the steady-state growth of a new 
phase in a strip geometry are obtained and different pattern formation processes in this system are 
investigated. 
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Many magnetic, superconducting and structural phase 
transitions in solids are accompanied by small lattice dis- 
tortions which lead to the presence of elastic deforma- 
tions. For each phase transition a characteristic defor- 
mation can be assigned, i.e., the deformation acquired 
by the new phase relative to the initial phase in the ab- 
sence of external forces. In some cases, these effects are of 
minor influence and can be ignored, but nevertheless for 
many applications the elastic strain causes qualitatively 
new and observable effects. The influence on the ther- 
modynamics of transitions between different phases has 
been thoroughly discussed in the literature (for a review 
see 0,0, and, for more recent developments, e.g. Q and 
references therein). 

One of the well known consequences is a thermody- 
namic elastic hysteresis, i.e. the splitting of the phase 
equilibrium point into two points, the points of the di- 
rect and inverse transition. It is mainly due to the co- 
herency at the interphase boundary, meaning that the 
lattice layers remains continuous through the boundary. 
Correspondingly, the hysteresis disappears without inter- 
face coherency |4j. Despite of the general understanding, 
some features of such transformations are still unclear, 
or at least under debate. For example, a distinctive two- 
phase equilibrium is established in the system within a 
certain temperature interval. The nature of this phe- 
nomenon is difficult to understand from the standpoint 
of ordinary thermodynamic equilibrium concepts: during 
the process of martensitic transformations the composi- 
tion of the phases does not change, and thus in such 
systems only one phase can be stable at a given temper- 
ature. 

However, the systematic theoretical study of the 
growth kinetics of such phase transitions accompanied 
by a lattice strain is much less advanced. In real systems 
the influence of elastic strain is often screened by many 
other effects, for example, by inhomogeneous composi- 
tions and temperature distributions, the Mullins-Sekerka 
instability, crystal anisotropy, polycrystalline structures, 
etc. Here, a phase-field modeling of such complicated 
systems can lead to qualitative descriptions of the kinet- 
ics of phase transitions in solids (see, for example, and 



references therein). The main purpose of this Letter is 
to develop a "minimum" kinetic model from which even 
quantitative results concerning the influence of strain ef- 
fects can be obtained. 

We start from the thermodynamical description of our 
model. The free energy density of an initial phase is 

F x =F* + ±ul + iml, (1) 

where F® is the free energy density without elastic effects, 
Uik is the strain tensor, A and /i are the elastic moduli 
of isotropic linear elasticity. The free energy density of a 
new phase is 

F2 =F% + \ («« - ul) 2 + n (u lk - ulf , (2) 

with u® k being a characteristic lattice strain assigned to 
the phase transition. 

Let us consider the simplest case, u® k = edik, at first. 
We assume that the elastic effects are small, e <C 1, and 
neglect the difference between the clastic coefficients in 
the two phases. Since in our description the reference 
state for both phases is the undeformed initial phase 
(see Eqs. and ©), the coherency condition reads 
uv 1 ) = u( 2 \ where u is the displacement vector. The 
superscripts (1) and (2) refer to the initial and the newly 
created phase respectively. Mechanical equilibrium at the 
interface demands o4n = c^nn and a4V = Out • Here, the 
stress tensor is given by a ik = ^(dF/duik + dF/duki); 
the indices n and r denote the normal and tangential di- 
rections with respect to the interface. The condition of 
phase equilibrium requires the continuity of a new poten- 
tial F across the flat interface [(J, 

F F ^nn^nn 2cr nr li nr . (3) 

In the general case of curved interfaces also the surface 
energy 7 should be taken into account, and the phase 
equilibrium condition reads in the case of isotropic sur- 
face energy 

F x - F 2 - 7 /C = 0, (4) 
where 1C is the curvature of the interface. 
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FIG. 1: Shapes of propagating fingers calculated for three 
values of the driving force. Top: A = 0.5 < At. Center: 
A = 0.85 > A t . Bottom: A = 1.05. The Poisson ratio is 
u = 1/3 and L/L* = 10. 



Eq. we find the volume fraction of the new phase, 



1 -2v 
l + v 



(7) 



which defines a dimensionless driving force A for this 
process. Then the parameter range for coexistence is 
< A < 1. The total energy gain of this two-phase 
configuration compared to the unstrained initial phase is 



AT = L[F? - cF 2 - (1 - c)fl] - 2 7 . 
Finally, using Eq. J7J we find 

A_ 2 

A? 



(8) 
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A critical nucleus of the new phase inside an un- 
bounded initial phase can exist only if 



(5) 



where E and v are Young's modulus and Poisson ratio 
respectively. This condition corresponds to the elastic 
hysteresis mentioned above. It can be obtained using the 
analogy of this elastic problem to the problem of thermal 
expansion for a given temperature field, as described in 


We discuss a simple strip configuration which allows 
the steady-state growth of an elastic "finger" consisting 
of the new phase (see Fig.nj. The unstrained elastic strip 
of width L is attached to fixed grips at the upper and 
lower boundary (u = there) and initially composed of 
the reference phase. We discuss a two-dimensional elas- 
tic problem using plane strain conditions (u z = 0). Also, 
we assume the complete wetting of the walls by the ini- 
tial phase. Thus, the new phase avoids a direct touching 
of the walls (see Fig. QJ. Far ahead of the propagating 
finger the initial phase remains unstrained. In contrast, 
far behind the tip a phase coexistence is possible within 
a certain parameter interval near the transition temper- 
ature, which is due to elastic effects. In this region the 
only nonvanishing component of the displacement vector 
is u y . The strain tensors are constant in both phases and 
their nonzero components are connected to each other by 



the relation (1 — c)u 
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= 0, in order to fulfill the 



conditions u y — on the walls and u 



(i) _ 
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= u y 2 ^ at the in- 
terface. Here, c is the volume fraction of the new phase. 
Then, the mechanical equilibrium condition gives 
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Taking into account the phase equilibrium condition 



where 



L* = 4 



(l-2i/)(l-i/) 7 
(1 + v)Ee 2 



(10) 



The finger grows if AT > 0, or, equivalently, if A > A^. 

Two remarks are in order. First, we obtain the two- 
phase structure only because we take into account the 
elastic effects and use the fixed volume boundary condi- 
tion. For stress free boundaries, we would obtain c = 1 
for any driving force F° — F° above the threshold (J5J and 
c = below the threshold. Second, the value of the driv- 
ing force for the transition can be controlled not only by 
temperature, but also by external strain. In particular, 
if the strip of width L is stretched by 6L, the homoge- 
neous term SL/ L has to be added to the strain u yy . This 
eventually leads only to a renormalization of the driving 
force, 



f?-F°-F°-i? 
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(11) 



Following our general aim to develop a minimum ki- 
netic model, we assume that the growth is controlled only 
by interface kinetics. Then, the local equation of motion 
of the interface reads 



k{F 1 -F 2 - 7 /C), 



(12) 



where v n is the normal velocity and k a kinetic coefficient. 

The conservation of energy requires that the excess 
AT is compensated by dissipation at the interface. This 
leads to a relation between the growth velocity and the 
driving force for the process, 



v 

VqL 



nlds - A 2 - A 2 , (13) 



where v is the steady-state velocity of the finger and 

k(1 + v)Ee 2 



2(1 - 2z/)(l - v) 



(14) 
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is the characteristic velocity scale for this system; n x is 
the projection of the interface normal on the growth di- 
rection x and the integration is performed along the in- 
terface. The dimensionless quantity J n~.ds/L is a com- 
plicated function of the parameters A, A^, v\ in the 
case of dynamical elasticity it also depends on the ratio 
vq/c s where c s is the shear wave speed. Near the equi- 
librium point A = Ai, the growth velocity behaves as 
v/vo oc A — Al. These results are valid for the dimen- 
sionless driving force A < 1. For A > 1 the fraction of 
the second phase becomes unity and Eq. I|13|l should be 
replaced by 

^ J n 2 x ds = 2A-1-A 2 L . (15) 

In order to obtain quantitative results for this problem 
we use a phase field code together with elastodynamics 
to describe phase transformations accompanied by stress, 
which we developed recently @. Let <fi denote the phase 
field with values <f> = 1 for the initial phase and = 
for the new phase. The energy density contribution is 
F = Fih(<f>) + i<2[l — h((f))], where the switching func- 
tion h{4>) — (f> 2 (3 — 2(f)) interpolates between the phases. 
The gradient energy is F s ((f) — 37£(V0) 2 /2 with the 
interface width £. Finally, F^ w = 670 2 (1 — 4>) 2 /£ is the 
double well potential. Thus the total energy functional is 
F = J dV (F + F s + Fdw)- The elastodynamic equations 
are derived from the energy by variation with respect to 
the displacements m, pili — —SF/Sui, where p is the 
mass density. The dissipative phase field dynamics fol- 
lows from <ft = — (n/3£)6F/5(j). These equations lead in 
the limit £ — > to the sharp interface description above. 
For the case of static elasticity, this was proved in ||. 

For the numerical realization, we employ explicit rep- 
resentations of both the elastodynamic equations and the 
phase field dynamics. The elastic displacements are de- 
fined on a staggered grid 0|- We shift the grid hori- 
zontally in order to keep the propagating tip always in 
the center of the strip; this allows to study steady-state 
growth in moderately large systems. The intrinsic length 
scale L* is chosen to be larger than the phase-field in- 
terface width, L* = 8£, and £ = 5 Ax, where Ax is the 
numerical lattice unit. We have performed calculations 
for the Poisson ratio v — 1/3 and for two values of the 
dimensionless strip width L/L* = 10 and L/L* = 20, or 
equivalently, for two values of the parameter A^ « 0.32 
and Ai « 0.22. All simulations are conducted on the 
parallel computer JUBL operated at the Research Cen- 
ter Jiilich. 

First of all, we have checked that the asymptotic 
conditions in the two-phase region far behind the tip, 
Eqs. 10 and J7J, are reproduced by our numerics with 
high precision. Though the appearance of the Asaro- 
Tiller-Grinfeld instability 0, [l3| might be naively ex- 
pected in the tail region because of the presence of non- 
hydrostatic stresses, this effect is never observed in the 
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FIG. 2: The dimensionless growth velocity v/vo as a function 
of the dimensionless driving force A. 

simulations. In the limit L — > oo, we have checked analyt- 
ically that the system is stable. We have also confirmed 
the predicted renormalization of the driving force l|ll|) 
due to an external strain SL/L. Next, we have investi- 
gated the properties of growing elastic fingers for different 
driving forces, F® — F®. After a transient regime the fin- 
ger always reaches the steady-state configuration. Char- 
acteristic shapes of the stationary fingers are presented in 
Fig-d The finger selects a symmetrical shape even if the 
initial configuration was vertically off-centered, i.e. the 
symmetrical configuration is stable. We note the ex- 
istence of the "bubble" in the tip region (top panel of 
Fig. [IJ. The bubble touches the walls at a specific value 
A = At < 1. With further increase of the driving force 
the touching region smoothly increases and diverges at 
A = 1 (see middle and bottom panel of Fig. ^1 . 

The dependence of the dimensionless steady-state 
growth velocity v/vq on the dimensionless driving force 
A is presented in Fig. [21 The velocity rises from zero 
at A = Al with increasing driving force. It turns out 
that the dependence is almost linear in a wide range, 
Al < A < Af. This is a slightly unexpected result, be- 
cause the linear behavior is anticipated to be valid only 
in the close vicinity of the equilibrium point A = A^. 
At both critical points, A = A t and A = 1, the growth 
velocity is a continuous function of the driving force. 

All the results presented so far have been obtained for 
the characteristic velocity scale vq being much smaller 
than shear wave speed, i.e. in the limit of static elasticity. 
We have performed additional runs with higher charac- 
teristic velocities such that v/c s ~ 1/2, where dynami- 
cal effects are relevant. For the same driving force, the 
dimensionless velocity v/vq decreases compared to the 
static elasticity limit, and the size of the bubble grows in 
order to reach the same dissipation according to Eq. I|13(l . 

Let us shortly discuss another simple example of tran- 
sitions involving shear strain in hexagonal crystals. For 
the transitions lowering the symmetry from Cq to C2 the 
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shear strain in a basic plane appears. For simplicity we 
neglect all other possible strains with higher (axial) sym- 
metry. We assume that the crystal is attached to two 
parallel walls as before. Let the principal axis Cq be 
oriented in z direction. By proper choice of the crystal 
orientation around the main axis in the initial phase, we 
obtain the new phase in three possible states having the 
following nonvanishing components of the strain tensor 
?;° • 



£ cos 28, 



e sin 28, 



(16) 



where the angle 8 = 0,±27r/3. Because the elasticity of 
hexagonal crystals is axisymmetric in harmonic approxi- 
mation and u® z = Ui Z = for the discussed problem, we 
can use Eqs. and J2J) for the energy densities of the 
two phases (see e.g. 01 )■ The moduli of the effective 
isotropic elasticity, A and /i, can be expressed in terms of 
the elastic constants of the original hexagonal crystal. 

A straightforward analysis of the stress state far be- 
hind the tip shows that among the possible configura- 
tions of new phases the energetically most favorable sce- 
nario are bicrystals, 8 = ±27r/3, as presented in Fig. [21 In 
the asymptotic tail region u^ty = and u^y = =Fe-\/3/2, 
where different signs correspond to different domains of 
the bicrystals. The distribution of the strain component 
u yy and the fraction of the new phase c can be readily 
found in the same way as before. For example, Eq. l(7)l 
should be replaced by 



= A = 



1 



1-2// 



4(1 ~ V2) (F? F°) 1 



E£ 2 



(17) 



The presence of the twin boundary with interfacial en- 
ergy 7b requires also a modification of the characteristic 
length and velocity scales, 



L* = 



i (l-^)(2 7 + 7b ) 
(1 - 2v)E£ 2 ' 



vo = k(2 7 + lb )/L* 



compared to Eqs. J7J) and l|14l) ; moreover, it leads to the 
existence of a triple junction in the tip region. 

For a numerical phase field study of these twin struc- 
tures we immediately take into account the symmetry of 
the appearing patterns and describe only either the up- 
per or lower half of the strip. At the symmetry plane 



the boundary conditions are u y = 0, a x 



0, and, for 



the specific case 71 « 7 considered here, dfi/dy = 0. 
This avoids more complicated multiphase descriptions 
which are in principle capable to describe the three dif- 
ferent phases. Although both patterns in Fig. [21 are en- 
ergetically equivalent far away from the tip, symmetry is 
broken by the choice of the propagation direction. For 
growth to the right, the orientation 8 — 2n/3 in the upper 
and 8 = — 27r/3 in the lower half, leading to propagation 
with a single tip (upper panel of Fig. [2J), is superior to 
the opposite case with repelling fingers (lower panel), as 
the growth velocity is higher, see Fig. 0] 




FIG. 3: Growth of bicrystal patterns for L/L* =5. The shad- 
ing illustrates the orientation of the shear strain u xy . Forward 
slashes correspond to 9 — —2n/3, backslashes to 9 — 2n/3. 
The strip lengths used in the simulations are much bigger 
than in the sections shown here; far away in the tail region 
both shapes have concentrations c — A — 0.6. The growth 
velocities are v/vq — 1.14 (top) and v/vq — 0.48 (bottom). 
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FIG. 4: Growth velocity of the different bicrystals for L/L* — 
5. A slight discrepancy from the equilibrium point Al is due 
to the accumulation of elastic energy in the transition region 
of the phase field, which leads to a renormalization of the 
surface energy. We checked numerically that this effect is 
suppressed in the sharp interface limit, as expected. 



In summary, a simple sharp-interface kinetic model of 
strain influenced phase transitions has been developed 
together with its phase-field realization. We obtained 
quantitative results for the steady-state growth of an elas- 
tic finger in a strip geometry and discussed the peculiar 
behavior of different pattern formation processes in this 
system. Influence of additional conserved fields, e.g. com- 
positional and temperature field, can be subject of future 
investigations. 
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